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We study the creation of a black hole (BH) pair separated by a domain wall, in the presence of 
a cosmological constant. We construct the solution representing a BH pair with a domain wall and 
compute the Euclidean action to evaluate the probability of the pair creation in the background with 
a preexisting domain wall. The BHs can be either neutral or magnetically charged ones. We compare 
the results of the charged case with those of the neutral case with the same cosmological constant. We 
find that the production rate of a charged BH pair is always suppressed in comparison with that of 
the neutral one in both four and five dimensions, irrespective of the sign of the cosmological constant. 
The Euclidean action is equal to the minus of the entropy. Since the horizon area of a BH is decreased 
as the magnitude of its charge is increased in general, the decreasing creation rate can be understood 
in terms of the increasing charge. We obtain the explicit confirmation on the relation between the 
pair creation rate of the charged BHs and the area of horizons in both the four- and five-dimensional 
cases in the presence of a cosmological constant. The singularity of the domain wall universe with 
charged BHs, as distinct from that with neutral BHs, can be avoided. 
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1 Introduction 



Black holes (BHs) are the fundamental objects in studying gravity, astrophysics and even high energy 
physics. The astrophysical BHs are the most promising sources of the gravitational waves and important 
for testing the general relativity in strong gravity regime. In cosmology, the primordial BHs will be probes 
for cosmic history and structure formations [1]. Furthermore, recent developments in higher-dimensional 
theories have brought our interests in finding new highcr-dimcnsional BH solutions as well as other kinds 
of black objects, and investigating the possibility of their productions in the high energy colliders as LHC 

[HISlll]. 

The vacuum in the presence of a strong field can decay through the so-called Schwingcr mechanism 
[5] , creating a pair of particles in an external electric field. Such particle production is one of the decaying 
processes of the background field in the given vacuum state (for recent works, see e.g., [SlIZ])- As the 
gravitational analogy of the mechanism, it was proposed that the pair creation of BHs, which represents 
a nonperturbative topological fluctuation of the gravitational field, could be possible in the background 
magnetic field [8]. The energy needed for the creation is provided by the background field. This was 
confirmed in [5] by constructing an appropriate Ernst instanton solution which describes a pair creation 
of BHs. The BH pair creation involving a cosmic string has been studied [TUl El HH US] ■ In this case, the 
energy to create the BH pair comes from the string tension. The pair creation has been also extensively 
investigated in the early universe [Mj [T5j [161 [17] , where the energy comes from the cosmological constant, 
and in other frameworks [IHl [HI [101 [HI [12]- In the case with a domain wall, the wall's gravitational 
repulsive energy [53] can give rise to a pair creation of BHs. Such a process was first studied in [53] using 
a cut-and-paste method. The topological BHs in the presence of the domain wall [55], and the creation 
of an inflationary brancworld with a pair of black cigars [26] have been studied. 

There are two kinds of approach to the creation of a BH pair. One is the bounce approach. In this 
approach, one needs to flnd an instanton solution of the Euclidean fleld equations, which interpolates the 
initial state without BHs and the flnal state with BHs. The other is the quantum cosmological approach. 
In this approach, one can employ the no-boundary proposal which represents the probability of the 
final state from nothing and that of the initial state from nothing. The creation rate from the initial to 
the flnal state can be obtained from these two probabilities. This method docs not need the interpolating 
instanton solution. In the semiclassical approximation, the formula of the production rate without a 
prcfactor in the bounce approach is equivalent to the formula in the quantum cosmological approach |28| . 

In this paper, we revisit the issue of the creation of a BH pair separated by a domain wall with 
an arbitrary cosmological constant. One of our motivations is the expectation that the creation of the 
domain wall with a charged BH pair in the flve-dimensional spacetime may provide a model of the 
braneworld universe. In terms of string theory or supergravity, the theory generically contains higher- 
rank antisymmetric tensor fields and BHs may be charged. Thus it is naturally expected that in any 
string-inspired model the domain wall universe is created together with the charged BHs. In addition, 
from the cosmological points of view, the possibility of the domain wall universe with charged BHs is worth 
being considered since the initially contracting universe can experience a bounce and the singularity can 
be avoided. From these points of view, it is very significant to evaluate the nucleation rate of the charged 
BHs. Thus, we could get insights on the quantum cosmological origin of braneworlds. We will construct 
the solution representing the charged and the neutral BH pair with a preexisting domain wall in the four- 
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and five-dimensional spacetimes. We will compute their production rates, and then compare the actions 
of the charged case with those of the neutral case with a same cosmological constant rather than obtaining 
the interpolating instanton solutions. 

The spacetime of the preexisting wall can be provided either by a cut-and-paste method or by the 
instanton solution mediating tunneling between the degenerate vacua in curved space [291 1301 1311 132j . The 
braneworld-like object can be obtained as the interpolating instanton solution by applying the mechanism 
in |31U32j . From this point of view, our present work can be the basic framework to make the spacetime, 
where we can describe the dynamics of a domain wall universe in the charged BH spacetime not only 
in the Einstein-Maxwell theory in Refs. [531 131] but also in the more general U{1) gauge theories in 
Refs. [35j|36l[32l- 

The organization of this paper is as follows: In Sec. 2, we review the magnetically charged BHs for 
general spacetime dimensions in the presence of a cosmological constant. In Sec. 3, we construct the 
solution for a charged BH pair separated with Z2 symmetry by a domain wall in the four- and five- 
dimensional spacetimes. The case of a neutral BH pair including non-Z2 symmetry was studied in j30) 
using the Bousso-Hawking normalization |15[ 116] in the bulk, where ambiguity of the time periodicity 
on the wall remains. In Sec. 4, we evaluate the Euclidean action and the resultant pair creation rate 
for the above system using the time periodicity by the Hawking temperature of BHs [38j in four- and 
five- dimensional spacetimes. We will show that the creation rate of a BH pair with the domain wall 
is proportional to the area of the horizons of created BHs in four- and five-dimensional spacetimes. We 
will obtain the first explicit confirmation on the relations between the pair creation rate and the area 
of horizons. And the application to the braneworld cosmology is discussed. In Sec. 5, we summarize 
our results and discuss the possible generalizations of our studies. It is natural to expect that the same 
properties are hold also in more than six-dimensional spacetime which will be discussed. 



2 Magnetic BHs 



We consider a system composed of Einstein gravity with either positive or negative cosmological constant 
A, coupled to the antisymmetric tensor field strength of rank (n — 1) F(^n-i) 
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where g = dctg^i^, q = dctqat, K and Kq are traces of the extrinsic curvatures of dM in the metric 
g^jy and 77^^, respectively, Fa^...a^_^ ~ [n — 1)V[q,j^^q,2...q^_ j , and is a unit normal vector which 
points outward to the boundary. The Greek indices run the (n -f- l)-dimensional spacetime, while the 
Roman indices run the 71-dimensional spacetime. The third term on the right-hand side corresponds to 
the boundary term for Einstein gravity [39]. The fourth term corresponds to the boundary term for the 
bulk form field. The term is not needed in the case of a magnetic field. In the electric case, there is the 
subtlety that the field must be purely imaginary on the Euclidean section and we must keep the boundary 
term for the bulk part [JD]. Note, however, the partition function in a definite charge sector for the electric 
cases with the boundary term is the same as that for the magnetic cases without the boundary term in 
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Figure 1: The phase of a magneticany charged BH with a positive cosmological constant in five- 
dimensional spacetime is shown. The vertical and horizontal axes show Aq and Am, respectively. The 
solid (red) curve corresponds to the extremal solutions = r_ (Eq. (O), while the dashed (blue) curve 
corresponds to the charged Nariai solutions r+ = Tc (Eq. ([5])). The intersecting point of these curves 
corresponds to the ultracold solution r+ = r_ = Tc- 



the semiclassical approximation. Hereafter we will omit the boundary term for the form field, because it 
will vanish in this paper. 

Varying the action with respect to the metric gives the Einstein equation 
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R,. = ^9,. + [in - l)(i^[Li]),, - ^.9m.(^^[Li])J , (2) 

where we have defined {F^^^_^^'^ = F^"^"'°'("-^'i Fi,ai---a^„_2) and {F[n_^) is its trace. 

A magnetically charged BH solution in the (anti-) de Sitter ((A)dS) spacetime is given by 

as =-f{r)dt + +7- c^^2(„_l^, /(r) = 1 - 
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where eai - a„_i = ±1, 7(n-i) is the volume element of the (n — l)-sphere, and indices {a^} are those for 
it. The other components of the form field are zero. The existence and number of the horizons crucially 
depend on the parameters. In this paper, we will focus on the cases of four- and five-dimensional BHs. 

2.1 The four-dimensional case 

In the case of n = 3, the solution Eq. ([3]) becomes the magnetic Reissner-Nordstroem-(A)dS BH in four 
dimensions 

ds' = -f{T)dt^ + ^+ r\de^ + sin' ed<p'), fir) = 1-^-^ + 4, 
jir) S r 

Fe0 = — sin0. (4) 



There are at most two BH horizons r± (r_|_ > r_). There is also a cosmological horizon outside the 
BH horizons Vc > ?>. The phase is quite similar to Fig. 1 for A > (see also [TB]) and Fig. 2 for A < 0, 
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Figure 2: The phase of a magnetically charged BH with a negative cosmological constant in five- 
dimensional spacetime is shown. The vertical and horizontal axes show \ls\q and |A|to. The solid curve 
corresponds to the extremal solutions r+ — r^. For solutions corresponding to the side below this curve, 
there are two horizons, and for those left to this curve there is no horizon. 



replacing Am and Ag with krr? and Aq^ respectively, where Figs. 1-2 represent the five-dimensional cases. 
The solid (red) curve in Fig. 1 corresponds to the extremal solutions of r_|_ = r_, while the dashed (blue) 
curve corresponds to the charged Nariai solutions of = re- The point where these two curves meet 
represents the ultracold solution of r_|- = r_ ~ Vc- In the region surrounded by these two curves and q = 
axis, three horizons exist at the same time. For the case of g = to, obtained from /'(r+) = — /'(rc), there 
is the special solution with three horizons, known as the lukewarm solution. The horizon positions in this 
solution are simply given by 

Such a simple analytic expression for the lukewarm solutions can not be found in other dimensions, 
although the lukewarm solutions themselves can also exist. 

For A < 0, the phase diagram is essentially the same as the Fig. 2, replacing |A|to and \A\q with |A|to^ 
and lAj^/'^, respectively. The solid curve shows the extremal solution and below this curve a BH has both 
the outer and inner horizons. 



2.2 The five-dimensional case 

In the case of n = 4. the Eq. becomes 

j{r) 6 

I? ^1 -2-0 /a\ 
i'xB'l}^ sm xsmt', (6) 

where dn\ — dx^ + sin^ x{dS^ + sin^ 9d(fP) denotes the unit 3-sphere and we take e^oip = +1- 

For A > 0, the existence of three horizons is ensured by three curves on the (Am, Aq) plane, namely 
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(7 = and 

qA = fetA ^((-2 + 3Am) + 2(1 - Am)i)^ (7) 

qA = ^cnA := ^((-2 + 3Am) - 2(1 - Am)^)^ (8) 
V3 

where gcxt and gcN correspond to the charges of the extremal and charged Nariai solutions, which are 
r_|_ = r_ and 7'c = r-|_ , respectively. The phase of a magnetically BH with a positive cosmological constant 
in five dimensions is shown in Fig. 1. goxt and (/cn are shown by the solid (red) and dashed (blue) curves, 
respectively. There is no essential difference from the four-dimcnsional case. The point where these two 
curves meet, where Am = 1 and Aq — represents the ultracold solution, where r^^ = = Tc = \J^- 
There is also the lukewarm condition satisfying f'{r^) = —f'{rc). As mentioned previously, there is 
no simple analytic expression for the horizons as in four dimensions. Note that the ultracold solution 
corresponds to the common end point of the previous two curves. In the closed region surrounded by 
these curves, there are three nondegenerating horizons. Outside this region, there is only the cosmological 
horizon. The essential picture remains the same in other dimensions. 
For A < 0, the condition for existing two BH horizons is given by 

g|A| < QextlAI := -^(2(1 + |A|m)3/2 - (2 + 3|A|m))'/'. (9) 

The phase of a magnetically charged BH in five dimensions is shown in Fig. 2. Note that there is no 
endpoint of the curve. The essential picture remains the same in other dimensions. 

3 Charged BH pair separated by a domain wall 

After giving BH solutions in the four and five dimensions, we construct the solution representing a pair 
of charged BHs separated by a domain wall. We assume the case with the Z2 symmetry with respect to 
the domain wall. 

3.1 BH pair separated by a domain wall 

From now on, we work in the Euclidean section. We consider a pair of charged BHs separated by a domain 
wall. We focus on the Euclidean action involving the contribution of a boundary domain wall 

1 /„ . .\ 1 
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^2 ' (71—2)!^ "'[ai^'-a2---a„-i\ > (10) 

where M and dM represent the bulk spacetime and the domain wall, respectively. Here, a denotes the 
tension of the domain wall. The last two terms denote the boundary term for gravity and the form field. 
The magnetically charged BH solutions in the Euclidean section are given by 

dsl^f{r)de + ^^^r-d^l_,,, /(r) = l-^|^-^ + ^, (11) 

where the avoidance of the conical singularity at the outer horizon r = r_|_ imposes the periodicity of t 
coordinate to be 
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Here, /3 is the inverse of Hawking temperature of BHs. Applying the cut-and-paste method or employing 
the instanton solutions mediating tunneling between the degenerate vacua, the solution representing an 
oppositely charged BH pair separated by a domain wall is constructed. Now wc will employ the Israel 
junction condition [41j. In this framework, Einstein equations should be solved on either side of the 
wall. In the Euclidean section, a domain wall is moving along the trajectory characterized by the affine 
parameter t, r ~ i'{t). t is normalized to satisfy + ft = 1, where "dot" means the derivative with 
respect to r. The equation of motion for r(r) is given by — V{r) = 0, where the effective potential is 
given by 

^^^^ 

By redefining the parameters, the potential is given by 

nr) = -Ar^ + l-^ + ^, (14) 



where 



2A AtV^ 

■~ n{n-l) 4(n- 1)2- ^ ' 



From Eq. (jlSp . in general we find 
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A>^^y (16) 

n[n — 1) 

On the other hand, for A < Q (and A < 0), 

lAl < -3^. (17) 

n[n — 1) 

We then regard A as an input parameter, rather than the tension cr itself. A stationary domain wall 
can exist at an equilibrium point r — which satisfies 

n..) = ^|_=0. (18) 

The corresponding mass of the BH is denoted by m*. In the next subsections, we will obtain the explicit 
expressions for and in the four- and five-dimensional spacetimes. The induced metric on domain 
wall metric is given by 

dsl^^dr'+rldnl_,. (19) 

For the domain wall at the equilibrium point, the Euclidean proper time on the wall has the periodicity of 
\fY(T*)P^ which can be interpreted as the inverse of the temperature measured by an observer living on 
the domain wall. The system composed of a domain wall together with a BH pair is now in the thermal 
equilibrium. 

In the semiclassical approximation, the production rate for a creation of a charged BH pair is given in 
Eq. (|10p. The boundary term for the form field vanishes for our solution. From now on, we rewrite the 
Euclidean action useful for the evaluation. By employing the Israel junction conditions [K] = —-p^K^a 
and the Einstein equations, the Euclidean action can be rewritten as 
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Note that the EucHdean geometry is smooth at the horizon because of our choice of /3 as Eq. (fT2)) . and 
hence there is no contribution of the horizon to Se- We therefore take only the boundary term of the 
domain wah into considerations. Our strategy is to compute the Euchdean action for our system in four- 
and five- dimensional spacetimes and compare our results of the charged case with those of the neutral 
case with a same cosmological constant. After plugging F^n^i-j = |;."~2(i-i) into Eq. (PH)) . which further 
reduces to 
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From Eq. ([TH]) and Eq. (HI]), wc find 
fir.) -{a- 

respectively. Thus, we finally obtain 
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where m* is the mass of the BH when the domain wall is placed at r = r*. For later convenience, we 
introduce the dimensionless quantities measured in the unit of charge 

-, a:=q^^c7. (24) 
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3.2 The four-dimensional case 



For the case of n = 3, from Eq. (|T8|) we obtain 



-l + Ari 
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For A > 0, the solution to Eq. (|25|) is given by 

1/2 

r*.± = I — ; I , m*,± 



0, 2Ari 
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V2 1 + 12i± Vl - 12i 



(25) 



(26) 



Thus, for < i < ^ 



6V3yl y^i±v/l-12i 

The positions of r*,_|_ and r*^_ correspond to the local maximum and minimum of V(r), respectively. 

^2 , there are two roots, which we call the (+)- and (-)-branch, respectively. For 

± > l?^- For a (+)-branch domain wall, + > 1^ 
for < A < For a (-)-branch domain wall, ^q^ < ml _ < for < A < yj- ^ special case \s A = 
where the degenerate root is given by = v^?, resulting in ml = ^q^ = |r^. This is the case of a pair 
of the ultracold BHs r_|- = r_ = Tc, where to^A ~ | and g^A = j. Thus, from the definition of A, we find 



a — Q. The other special case is a pair of the lukewarm BHs Eq. ([S]). This is obtained for 
namely, for A ^ -j^ of the (-l-)-branch and for A = of the (-)-branch. In the limit of g — > 0, 
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Figure 3: The mass of the BH m* is shown as the function of \A\ in the five-dimensional case. The sohd, 
dashed and thick (red) curves represent m*,+ , m*,_ (for A > 0) and m* (for A < 0) as functions of \A\, 
respectively. For the case of the four dimensions, we obtain essentially the same picture. There is no 
endpoint for the curve of A < 0. 



Thus, only for the (+)-branch can the domain wall exist. The corresponding mass of the BH is given by 

TO* I ^ — ]=. (28) 



Until now, we have not taken the existence of the outer horizon into consideration. However, it is 
straightforward to confirm that there is no real solution of the outer horizon for the (-l-)-branch of A < ^ 
and for the (— )-branch of A > It means that for the (— )-branch no horizon can be formed for A > ^ 
obtained from Eq. (jl6p with rt = 3. For the special lukewarm solution, r*,_|- = 2q = 2?7i*, and then 
Tc > > ?'-(_, since > 2q and r-|_ < 2q, which ensures that the (+)-branch domain wall is always located 
between the cosmological and outer horizons. On the other hand, however, the (— )-branch domain wall 
for the lukewarm solution may exist for A ^ with q ^ 0. But it contradicts the condition A > and 
hence no (— )-branch BH pair is formed. 

For A <0, the solution to Eq. (^5]) is given by 

1/2 



l + 12\A\\ ^ ^ -1 + 12|A| + ^1 + 12|A| 

^ — • (29) 



1 + Vl + 12|^| 



The position of corresponds to the local minimum of V{r). For \A\ -> 0, to* — > q and for larger \A\, 
TO* is increasing. In the limit of g — > 0, ?'* — > 0. 

However, it is straightforward to confirm that for A < no horizon is formed for \A\ < But from 
Eq. ([T7)) with 71 = 3, we must have |A| < -ly-. Therefore, the horizon cannot be formed for A < 0. 

The phase diagram for a domain wall in four dimensions is quite similar to Fig. 3 which is for the 
five-dimensional case discussed in the next subsection. 
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3.3 The five-dimensional case 



Along the similar arguments to the case of four dimensions, we discuss the properties of a domain wall in 
the five-dimensional spacetime. For A > 0, the position of the domain is given by either of the following 
two branches 



1 r 
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The positions of r*,_|_ and r*__ correspond to the local maximum and minimum of V{r), respectively. For 
A < 0, 



f2 



1 



6\A\ 



-l + D{\A\y 



-1/3 



D{\A\) 



1/3 



(31) 



Here, the position of corresponds to the local minimum of V{r). We defined 

C(i) 1 - 54i2 + 6\/-3i2 + 8li4, 



D{\A\) -1 + 5A\A\' + 6^J -i\A\^ + 81\A\\ 



(32) 



and the expression for r*,_ is real. Thus, we find A < j^- In the limit of g — > for A > 0, only for the 
(+)-branch there is the domain wall solution ^ 



jj. For A < 0, there is no domain wall solution of 



-> 0. 

For A > 0, the corresponding BH mass is given by 



m^^+{A) 
m»,_(A) 



^(1 - (7-2/3(i) + 2(7-i/3(i) + 2C^/^{A) ~ C2/3(i)), 
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2[C(i)i/3(3(l - i%/3)V-3i2 + 8li4 + C(i)i/3) 
27i2(2(l + iV3) - C(i)i/3(i _ iV3))] |. 



Similarly for A < 0, the mass becomes 



m,(i) = - 1 + D-^^^{A) + 2D-i/3(i) + 2L>i/3(i) + D^/^{A)). 



(33) 



(34) 



In Fig. 3, we have shown the phase of a domain wall in five dimensions. For A ~ I^^j mass in the both 
branch becomes the degenerate value to* = ^q, which corresponds to the ultracold BHs. Then, from 
Eq. (fTS)) for n = 4 with A = of the ultracold case, it turns out that the corresponding domain wall is 
tensionlcss. In the limit of q — > 0, the mass of the domain wall for the (+)-branch becomes to*.-|_ — > 

As in the four-dimensional case, however, it is straightforward to confirm that there is no outer horizon 
for the (+)-branch of A < ^ and of for the (— )-branch of A > ^. Therefore, for the (— )-branch no horizon 
can be formed since A> ^ from Eq. (fT6)) with n = 4. Similarly, it is straightforward to confirm that for 
A < no horizon can exist for |^| < From Eq. ([TT]) with 71 = 4, we must have \A\ < ^. Therefore, 
there is no BH pair for A < Q. 
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4 Pair creation rates 



In the semiclassical approximation, the formula of the production rate without a prefactor in the bounce 
approach is equivalent to the formula in the quantum cosmological approach |28j . In the bounce approach, 
the production rate for a neutral BH pair in the background spacetime with a domain wall is given by 
Ti = j\(.-lSEinbis)~SE{bg)]^ where Ssinbls) and Ssibg) mean the Euclidean action of neutral BHs with a 
domain wall and the action of the wall, respectively. The charged case also has the same form. Thus the 
ratio of two production rates is written by F = T2/T1 = p-lSE{cbis)-SE{nbis)] ^ where SEicbls) means the 
Euclidean action of charged BHs with a domain wall. In the quantum cosmological approach, the ratio of 
two probabilities is given by T = P2/P1 = e-[SE(cf'is)-SB(nWs)] ^ .^^j^j,^.^ ^ ^-SEicbis) ^ ^-SE(nbis) 

mean the probability of charged BHs with a domain wall from nothing and of neutral BHs with a domain 
wall from nothing. This method does not need the interpolating instanton solution. In the approximation, 
the half of the Euclidean action is used in the wave function. The probability is related to the real part 
of the Euclidean action in the quantum cosmological approach. This action is equivalent to the action, 
Se, in the present work. 

In this section, wc compute the Euclidean action for each solution representing a BH pair separated by 
a domain wall in four- and five-dimensional spacetimes by employing the quantum cosmological approach. 
As wc have seen in the previous section, the outer horizon of BHs can be formed only for the (+)-branch 
of the A > case. Thus, in this section we will focus on these cases. We then evaluate the nucleation 
rates. 



4.1 The four-dimensional case 
4.1.1 For a pair of neutral BHs 

For a pair of the four-dimensional neutral BHs, a domain wall solution can exist only for A > 0, with 
m^, = ^J^^ and = 3m* = "^/l^- Thus, we have to focus on the case of A > 0. From Eq. ([25]) for n = 3 
with we obtain 

From now on, we discuss the cases of A > and of A < 0, separately. For A > 0, we obtain 

r+ = ^cos(i(7r + 0)), (36) 

where cos9 := and sm9 = ^/l — x with x := = Here a quantity with "tilde" is dimensionless, 
and normalized by an appropriate power of A. Since < a; < 1, < cos (i(7r -t- 6*)) < i. Finally, the 
Euclidean action is obtained by 

^F4,5ch(i), Fi,sch{A) := 3cos2 (i( 
To derive it, we have employed the formula cos3q! ~ 4cos'^ a — 3 cos a. It is easy to confirm that 



SE{nbls) = ^-—Fi,sch{A), F4,5c/«(^) :=3cos2(-(7r + 0)). (37) 



which leads to 



3r2_A 

FiSch — — -. — , (38) 



SEinbls) ^ -—^ ^ ^ --, (39) 
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where A ~ 2 x (47rr^) and G = ^ represent the area of horizons (the factor 2 means two sides) and the 
efifective gravitational coupling. Thus, the absolute value of the Euclidean action is equal to the entropy 
of created BHs, following the area law. This result is somewhat natural since the domain wall with a pair 
of BHs is now in the thermal ec^uilibrium. As we will see below, the same relation also holds for a charged 
BH pair with a domain wall, as well as in the five-dimensional cases. Note that the cosmological horizon 
is always out of our system and the thermal contribution from it is absent. 
Similarly, for A < we obtain 

where sinhG := with U = jj- Here a quantity with tilde is dimensionless, and normalized by 

an appropriate power of |A|. Finally, the Euclidean action is obtained by 

SEinbls) ^ -^^^Ti.schiA), J-4,sch(i) :=3sinh2(ie). (41) 

To derive it, we employed the formula sinh(3Q!) = 4sinh^ a + 3 sinha. It is easy to confirm that 

F - (42) 
■Ti^Sch- — , (42) 

which leads to Eq. ([59]) . The absolute value of the Euclidean action is equal to the entropy of created 
BHs, following the area law. 

4.1.2 For a pair of lukewarm BHs 



For a charged BH pair, the only analytically tractable case is the case of a pair of lukewarm BHs in four 
dimensioi 
given by 



dimensions, in which ~ 2to* ~ ^^73' ^^"-^ satisfied. In this solution, the position of the outer horizon is 



From Eq. (gS]) for n = 3 with 



, ' .. (44) 



3A 



the Euclidean action becomes 



SEicbls) = -^F4,„.e(^), F4,..e(^):-Y^(l-^l-(^)'n'. (45) 



For a given A, it is straightforward to see F4^iuke{A) < F4,sch{A) for any A (See Fig. 4), where F^^sch 

is given in Eq. ([57)) . Thus, the production probability P2 = e"'^^^'^'''*) of a lukewarm BH pair is always 

smaller than that of the neutral case. It is easy to find 

3rlA , , 

FiMke = (46) 

which leads to 

SEicUs) = -—^ = -^ = --. (47) 

As for the neutral case, the absolute value of the Euclidean action is equal to the entropy of created BHs, 
following the area law. 
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Figure 4: Fi{A, q), F^^^schiA) and F^^iukeiA) arc shown by the sohd (blue), dashed (red) and thick (green) 
curves, respectively, as the function of A for q = w 0.289. 



4.1.3 For a pair of general charged BHs 

We focus on the case of the (+)-branch domain wall with A > 0. Substituting rt = 3 into Eq. (123 



SE{cbls) 



167r^r4 



1 



2q2 

7^ 



1 



-r| A + rlA + q^[ 



1 



1 



(48) 



For A > 0, for comparison with the case of a neutral BH pair, it is useful to rewrite the Euclidean action 

(49) 



bE[cbls) = T . 



where 



FiiA,q) := 



3r+ 



1 



2f 



1 



3 



q ( — - — 



(50) 



Here a quantity with tilde is dimensionless, and normalized by an appropriate power of A. Note that 
F4(i,0) = F4.sch{A) defined in Eq. and FiiA, = F^jukeiA.) defined in Eq. (gS]). We then 



4^1/2, 

numerically evaluate the function F4{A,q) for the general cases. In Fig. 4, the comparison of F4{A,q), 
Fi^SchiA) and i^4,/„fee(^) is shown. Note that for this choice of parameters, i « 0.333 < A < 1. In Fig. 
5, F4 is shown as a function of q for a fixed A. We find that F4(A,q) < F4^sch(A). Therefore, the pair 
production rate for a charged BH is always smaller than that of the neutral case. We then numerically 
confirmed 

3rlA 



F4 



(51) 



which leads to Eq. (|T7)) . As for the neutral case, the absolute value of the Euclidean action is equal to the 
entropy of created BHs, following the area law. Then, the decreasing pair creation rate for the charged 
case can be understood in terms of the decreasing area of the outer horizon for increasing charge with a 
fixed cosmological constant. 
Similarly for A < 0, 

647r2 



SEicbls) 



-F4, 



(52) 
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Fa 

0.5 r 




0.0 0.1 0.2 0.3 0.4 0.5 

Figure 5: F4^{A, q) is shown by the sohd (red) curve as the function of q for A — 0.5. 



where 



3r4 



3m » 



2q^ 



Afl 



3 V r+ 



(53) 



Here a quantity with tilde is dimensionless, and normaUzed by an appropriate power of |A|. Note that 
Fi{A,Q) = F4^^sch{A) defined in Eq. (|4T|) . In Fig. 6, the comparison of ^4(^,9) and Ti^sch{A) is shown. 
In Fig. 7, is shown as a function of q for a fixed A. Therefore, the pair production rate for a charged 
BH is smaller than that of the neutral case. We then numerically confirmed 



-^4 = 



irl\K\ 



(54) 



which leads to Eq. (|T7|) . The absolute value of the Euclidean action is equal to the entropy of created BHs, 
following the area law. Similarly, the decreasing pair creation rate for the charged case can be understood 
in terms of the decreasing area of the outer horizon for increasing charge with a fixed cosmological constant. 



4.2 The five-dimensional case 

In this subsection, we repeat the similar computations for the five-dimensional spacetime. 
4.2.1 For a pair of neutral BHs 

Similarly to the case of four dimensions, in the case of the neutral BHs (7 = 0, the domain wall solution 
exists only for A > 0, with = and = -^^= ~ 2m* . For A > and A < 0, the BH horizon is 
located at 




respectively. From Eq. (j23p for n = 4, for A > we obtain 

('2\/3'l'^7r^ ~ / I 1~N 3/2 

SEinhls) = -^-^J^F5,sch, F5,5ch(A):= (l-^l-^J , (56) 
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Figure 6: J^4{A,q) and J-4,sch{A) are shown by the soHd (blue) and dashed (red) curves, respectively, as 
the function of ^ for g = ^ w 0.289. 
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Figure 7: J-a^A, q) is shown by the solid (red) curve as the function of q ior A — 0.5. 
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where A 



^, while for A < we obtain 



where A := j^. It is easy to confirm that 



Ssinbls) = 



1 \3/2 

^,,Sck{A) :=(^/l + _-l_ 



(57) 



(58) 



irrespective of the sign of A, where A = 2x ) is the area of black horizon. As for the four-dimensional 

case, the absolute value of the Euclidean action is equal to the entropy of created BHs, following the area 
law. The result agrees with that for a bulk BH pair of Ref. [26] . 



4.2.2 For a pair of general charged BHs 

We focus on the case of the (+)-branch domain wall with A > Q. Substituting n — A into Eq. 
Euclidean action reduces to 



SE{cbls) 



^2 /'4m, 



1 3 



'irlA + 6q' 



For A > 0, 



where 



9v^ 4m. 



SEicbls) = -F5, 

K^As 



3Aff - 



3p 



- 1 



6q^ 



, the 



(59) 



(60) 



(61) 



Here a quantity with tilde is dimcnsionlcss, and normalized by an appropriate power of A. Note that 
F5{A,0) ~ F^^schiA) defined in Eq. ([56| . In Fig. 8, the comparison of F5{A,q) and F^^schiA) is shown. 
In Fig. 9, F5 is shown as a function of q for a fixed A. Wc find that F5(A, q) < F^^schiA). Therefore, as 
in the four-dimcnsional case, the pair production rate for a charged BH is always smaller than that of the 
neutral case. We then numerically confirmed 



F 



which leads to 



SE{cbls) 



r^j_A2 

WW' 



G 



A_ 

4G' 



(62) 



(63) 



where A = 2 x (27r^r'|_) and G = ^ represent the area of horizons (the factor 2 means two sides). The 
absolute value of the Euclidean action is equal to the entropy of created BHs, following the area law. 
Then, the decreasing pair creation rate for the charged case can be understood in terms of the decreasing 
area of the outer horizon for the increasing charge with a fixed cosmological constant. 
For A < 0, we obtain 



SEicbls) 



(2V37r)= 
" /t2|A|i 



(64) 
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Figure 8: Fr^{A, q), and F^^schiA) are shown by the soHd (blue) and dashed (red) curves, respectively, as 
the function of A for q = 0.3. 
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Figure 9: Fq{A, q) is shown by the solid (red) curve as the function of q for A = 0.5. 
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Figure 10: J^5{A,q) and J-5.sch{A) arc shown by the sohd (blue) and dashed (red) curves, respectively, 
as the function of ^ for g = w 0.192. 



where 



MA,q) := [3if4^!± + 6g2f )1. (65) 



+ 



Here a quantity with tilde is dimensionless, and normalized by an appropriate power of |A|. Note that 
J^5(A, 0) = J^5,sch{A-) defined in Eq. ([57| . In Fig. 10, the comparison of J^4{A, q) and J^5,sch{A) is shown. 
In Fig. 11, is shown as a function of q for a fixed A. We find that J^5{A,q) < J-5.sch{A). Therefore, 
the production rate for a charged BH pair is smaller than that of the neutral case. We then numerically 
confirmed 

(V3)3 ^ ' 

which leads to Eq. (|63p . As for the neutral case, the absolute value of the Euclidean action is equal to 
the entropy of created BHs, following the area law. Similarly, the decreasing pair creation rate for the 
charged case can be understood in terms of the decreasing area of the outer horizon for increasing charge 
with a fixed cosmological constant. 

4.3 Application to the braneworld cosmology 

After the nucleation, the domain wall may evolve in the radial direction of the bulk spacetime. The 
particular case of interest is that of the five-dimensional spacetime, where the moving domain wall may 
be interpreted as our braneworld universe. Thus in this subsection, we will focus on the five-dimensional 
case. The cosmological equations on the domain wall have been studied in [251 

The behavior after the nucleation crucially depends on the sign of A and the branch of the domain 
wall. For A > 0, there are two branches of r*__|- and r*._ (see Sec. 3.3). They correspond to the local 
maximum and minimum of the potential V{r), respectively. However, as we have discussed previously, 
for the (— )-branch domain wall the BH pair cannot be formed for A > ^ obtained from Eq. with 
n — A. Thus, wc focus on the (+)-branch domain wall. In the (-l-)-branch, the domain wall is initially the 
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Figure 11: J^5{A,q) is shown by the sohd (red) curve as the function of q for A ^ 0.5. 



Einstein static universe just after the nucleation, but may evolve once the wah position deviates from the 
equihbrium point due to the small perturbations. Assuming that the domain wall evolves as r = a{T), 
where a(r) plays the role of the scale factor of the Robertson- Walker universe on the domain wall and r 
is now the proper time in the Lorentzian signature, the cosmological equation is given by + V{a) = 0, 
hence 



1 1 , 2m, 
^ = :tA + — r 



3 



(67) 



where X :— 3A becomes the effective cosmological constant [IHl EH- The second and third terms on 
the right-hand side are induced due to the mass and charge of the created BHs, respectively. The mass 
term behaves as the radiation in the universe [5^1 IH] , while the charge term behaves as the stiff matter 
with a negative energy density. Cosmological solutions of the domain wall universe in the charged BH 
background have been studied in e.g., Ref. |33| . There are two possibilities, namely the expanding 
(a > 7'*,+) or contracting (a < 7'*,+) domain wall. The first case corresponds to the expanding universe, 
approaching the de Sitter inflation with A, since the contributions of the mass and charge terms are diluted. 
The second case corresponds to the initially collapsing universe. However, in this case the domain wall 
does not run into the singularity and experiences a bounce, since there is the barrier in V{a) because of 
the charge q. 

For A < 0, the BH pair cannot be formed for \A\ < ^ obtained from Eq. (flT)) with n = 4, and we do 
not consider this case. 



5 Summary and discussions 



In this paper, we have studied the creation of a BH pair separated by a domain wall in the four- or 
five-dimensional spacetimes with a cosmological constant. We have constructed the solution representing 
a BH pair with a domain wall. In any dimensions, the solution which involves the domain wall outside 
the outer BH horizon can be formed only at the local maximum of the potential, which has a smooth 
limit to the neutral case. The singularity of the domain wall universe with charged BHs, as distinct from 
that with neutral BHs, can be avoided. 
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There are two kinds of approach to the creation of a BH pair. One is the bounce approach. The other 
is the quantum cosmological approach. In the semiclassical approximation, the formula of the production 
rate without a prcfactor in the bounce approach is equivalent to the formula in the quantum cosmological 
approach fIE\ . 

We have computed the Euclidean action for the above system in four- and five- dimensional spacetimes 
by employing the quantum cosmological approach. We then have compared our results of the charged 
case with those of the neutral case with a same cosmological constant. We find that the production rate 
of a charged BH pair is always smaller than that in the case of a neutral BH pair in both dimensions. We 
also have confirmed that the Euclidean action is always equal to the minus entropy of the created BHs, 
following the area law. We have given the first explicit proofs on the relation between the pair creation rate 
and the area of horizons of the created BHs in the four- and five-dimensional spacetimes. The spacetime of 
the domain wall can be provided either by a cut-and-paste method or by the instanton solution mediating 
tunneling between the degenerate vacua in curved space [221 EOl IHIl 1221 ■ The braneworldlike object can 
be obtained as the interpolating instanton solution by applying the mechanism in |31l I32j . From this 
point of view, our present work can be the basic framework to make the spacetime, where we can describe 
the dynamics of a domain wall universe in the charged BH spacetime not only in the Einstein-Maxwell 
theory in Refs. [221 [Mj but also in the more general U{1) gauge theories in Refs. [23 [23 [2Zj • With this 
motivation, the application to the braneworld cosmology has also been discussed. The four-dimensional 
geometry just after the nucleation is the Einstein static universe. The branch which is smoothly connected 
to the neutral case corresponds to the local maximum of the potential and eventually evolves, while the 
others are the local minimum. 

In this paper, we have employed the time periodicity by the Hawking temperature of BHs |38] so 
that no conical singularity appears at the horizon. The Euclidean geometry is smooth at the horizon, 
and therefore we only employed the boundary term by the presence of the domain wall. We omitted the 
boundary term for the form field, because it vanished in this paper. 

We briefly discuss the possible generalization and the results expected from the previous results. 

The first straightforward generalization is to the case of more than six-dimensional spacetimes. Even 
in such a case, the phase diagrams of the BHs and domain walls are very similar to Figs. 1-2 and Fig. 3, 
respectively. Similarly, the instanton for a pair of BH with a domain wall can be constructed. Note that 
the horizon inside the domain wall can be formed only at the local maximum of the potential V(r) for A 
satisfying Eq. p^ . Then, the creation probability is also similarly evaluated, which gives the essentially 
the same results as in the previous cases. In particular, the Euclidean action for the charged BH pair with 
a domain wall is expected to be 

SEicbls) = (68) 

where A ^ 2x r2„_ir"~"'^ is the area of the outer horizon, where the factor 2 means two sides. Following the 
area law, {~Se) coincides with the BH entropy, which is expected from our construction of the instanton. 
The nucleation rate e""^^ for a charged BH pair is always suppressed in comparison with that of a neutral 
pair for a given cosmological constant, since the area of the outer horizon decreases as the charge increases. 

Although it may be most plausible that two identical BHs are produced across the wall, it would 
be interesting to discuss the possibility without Z2-symmetry. This is the case where the cosmological 
constant in both sides is different and the domain wall is interpolating two different vacua. In this case. 
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in terms of the continuity of the magnetic field strength across the wall, two BHs must have the same 
amount of charge. In the case with Z2 symmetry, for a given cosmological constant A, charge q and 
domain wall tension a (or A), we could determine the positions of the outer event horizon r+ and that of 
the domain wall r^, and the BH mass m*, through the horizon condition = and two stationary 

conditions for the wall Eq. psp . We can extend this argument to the case without Z2 symmetry. For the 
given bulk cosmological constant A/, charge q, and tension of the domain wall a, where / specifies the 
side with respect to the wall, totally five quantities, i.e., the masses m/ and horizon positions r^j of BHs 
(four quantities) and the domain wall position must be specified by the five independent conditions. In 
our case, they are the horizon condition in each side //(r+./) = (hence two conditions), two stationary 
conditions as Eq. ((T8)) . and the continuity of f{r*) across the wall. The last condition is due to the 
requirement that the temperature measured by an observer on the domain wall must be unique. Thus, in 
contrast to the case with Z2 symmetry, even if a solution for the above set of equations exists, the masses 
of BHs are gencrically different. Otherwise, the system becomes overdetermined. It would be interesting 
to look for instanton solutions satisfying the above five relations, and evaluate the nucleation probability, 
although we now leave these issues for future studies. 
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